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ABSTRACT 

We present a new analysis of the long-period variables in the Large Magel- 
lanic Cloud (LMC) from the MACHO Variable Star Catalog. Three-quarters of 
our sample of evolved, variable stars have periodic light curves. We characterize 
the stars in our sample using the multiple periods found in their frequency spec- 
tra. Additionally, we use single-epoch Two Micron All Sky Survey measurements 
to construct the average infrared light curves for different groups of these stars. 
Comparison with evolutionary models shows that stars on the red giant branch 
(RGB) or the early asymptotic giant branch (AGB) often show non-periodic vari- 
ability, but begin to pulsate with periods on the two shortest period-luminosity 
sequences (3 & 4) when they brighten to K s » 13. The stars on the thermally 
pulsing AGB are more likely to pulsate with longer periods that lie on the next 
two P-L sequences (1 & 2), including the sequence associated with the Miras in 
the LMC. The Petersen diagram and its variants show that multi-periodic stars 
on each pair of these sequences (3 & 4, and 1 & 2) typically pulsate with periods 
associated only with that pair. The periods in these multi-periodic stars become 
longer and stronger as the star evolves. We further constrain the mechanism be- 
hind the long secondary periods (LSPs) seen in half of our sample, and find that 
there is a close match between the luminosity functions of the LSP stars and all of 
the stars in our sample, and that these star's pulsation amplitudes are relatively 
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wavelength independent. Although this is characteristic of stellar multiplicity, 
the large number of these variables is problematic for that explanation. 

Subject headings: galaxies: individual (LMC) — stars: AGB and post-AGB - 
stars: variables: other 



Introduction 



Stellar pulsation of giant stars appears to be a ubiquitous and important phenomenon— 
RR Lyrae and Cepheid variables form the basis for the distance scales we use. Miras and 
other long-period variables (LPVs) however, are not as well understood, largely because their 
cool and tenuous atmospheres are dynamic environments with a great diversity of molecular 
species forming and disassociating as the star pulsates. In recent years however, these stars 
have attracted increase d attention as m i cro-le nsing surveys of the Large and Small Magel- 



l anic C louds (OGLE — iPaczvnski et al. 



1994 ): OGLE II — llJdalski. Kubiak. & Szvmanski 



jl997h : MACHO — klcock et alJ jl997h ) have produced large catalogs of LPVs. Well- 



sampled light curves and excellent photometry give us an opportunity to better understand 
both the mechanisms behind long-period variables and the physical processes at work in the 
latest stages of stellar evolution. 

Before the wealth of data from micro-lensing surveys, LPVs were traditionally classified 
by the amplitude and stability of their y ariability in the V band (e.g. The General Catalog 



of Variable Stars — iKholopov et al.l (119961 )). In this scheme, stars with well-defined pulsation 



are classified as Miras if the amplitude of their variability exceeds 2.5 magnitudes in V, and 
as Semi-Regular Type a (SRa) stars if not. Those with multiple periods, unstable periodicity, 
or poorly expressed periodicity, are classified as SRb stars. The MACHO Project's survey 
of the Large Mag ellanic Cloud (LM C) revealed five parallel sequences of LPVs in period- 
luminosity space (jCook et al.lll996l ). pro mpting a clas s ificati on scheme that uses the period 
of pulsation as its primary discriminator. I Wood et al.l (119991 ) identified the cause of the first 
three period-luminosity sequences (denoted A, B, and C) as pulsation, and suggested that 
the two longest period sequences (E and D) could be attributed to binary systems. The 
stars in Sequence E showed the characteristic light curves of contact binary systems, and 
Sequence D stars — those with the longest periods — simultaneously exhibited at least one 
shorter period that was coincident with Seque nce B . This is likely to be the LMC equivalent 
of the "long secondary periods" described by iHouki (119631 ) for Galactic LPVs, although the 
periods that comprise Sequence D are on average three times shorter than the LSPs listed in 
Houkl (119631 ). IWood et al.l (119991 ) proposed that these stars are composed of accreting binary 
systems, with the long period caused by partial eclipses due to an unseen, dust-enshrouded 
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compamon. 



In this work t he LMC period-luminosity sequences will be named in the manner of 
Fraser et al.l (120051 ). f rom shortest to lon gest period: 4, 3, 2, 1, E, and D. We retain the 
names D and E from IWood et al.l (119991 ) , but rename his Sequence C to Sequence 1 and 
count up toward the shorter periods. This approach provides for a graceful way to accom- 
modate additional short-period sequences. Indeed, the use of 2MASS K s magnitudes as the 
luminosity indic a tor ca used a split in the original Sequence B — producing Sequence s 2 and 3; 
Kiss fc Bedding! (120031 ). A fifth sequence was identified by lSoszynski et al.l (J2004aj) through 
the examination of all significant frequencies of these stars instead of just the strongest 
frequency. 

In general, stars brighten and redden as they evolve along the Red Giant Branch (RGB) 
and the Asymptotic Giant Branch (AGB). Since the typical J — K s color of the stars in each 
sequence reddens as we progress from the short-period Sequence 4 to the longer period Se- 
quence 1, this suggests that evolution proceeds from shorter periods toward longer periods, 
at least for Sequences 1-4 (IFraser et al.l 120051 ) . In fact, the low lumino s ity bases of Se- 
quen c es 2, 3, and 4 are he avily populated by RGB stars (llta et al.l 120021 ; iKiss fc Bedding 
20031 . 12004! ; llta et al.llioO^ . Above the tip of the RGB, models of AGB stars that include 
the effects of mass loss co nfirm that stars continue their evolution to higher luminosities 
( IVassiliadis fc Wood! Il993l ). In period-luminosity space, the Mira and SRa stars are not 
clearly separated, with SRa stars found throughout Sequences 1-4. 



Soszynski et al.l (j2004al ) identified a more useful division for LPVs in the LMC and 



Small Magellanic Cloud (SMC) than the Mira/SRa/SRb system. In this system an LPV 
is classified as an OSARG (Ogle Small Amplitude Red Giant) if one of its three strongest 
periods falls onto Sequence 4, and as an SRV/Mira if not. This division separates LPVs 
into two groups with a variety of distinct properties, and it also shows that Sequence D is 
composed of two different populations: a dimmer population that covers a relatively broad 
period range, and a more luminous, redder population that sh ows a tighter period - luminosity 
relationship (the color changes described here can be seen in Fraser et al. ( 20051 )). 

Stars in Sequence E show "ellipsoidal" light curves JSoszvnski et al.ll2004bh . where the 
brightness modulation is due to the gravitational distortion of one member of a close binary 
system. These light curves exhibit dual minima of unequal depths, but the effect is small 
enough that most methods (including ours) find a period for these systems that is half of 
the orbital period. We refer to this period as the Fourier period and use it to distinguish 
Sequences E and D in our plots. When stars in Sequence E a re plot ted on the period- 



l umin osity diagram at their orbital period (as in ISoszynski et al.l (l2004bl ) and iDerekas et al. 



( 120061 )) they smoothly join Sequence D. This, along with a recent analysis of OGLE light 
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( I Adams et al. 



2006|). 



curves in ISoszvriskil (20071), sup ports the binary companion explanation for Sequence D 
put forth by Wood et al.l ( 19991 ). as does a radial velocity study of several of these stars 



In lFraser et al.l (120051 ) (hereafter Paper I) we used the MACHO and 2MASS magnitudes 
and colors to characterize LMC stars in each period-luminosity sequence, and classified the 
stars in each sequence as Miras, SRa, or SRb. At that time, only 52 percent of the stars 
in the color- and magnitude-defined sample had a well-determined period in the MACHO 
Variable Star Catalog. In this work, we expand the successfully analyzed stars to 93 percent 
of our sample, as well as consider their multi-periodic properties. We also use our results to 
describe the characteristic variability at each of the stages of RGB and AGB evolution by 
comparison with population synthesis models, including the stars that show very weak or 
non-existent periodicity. 



Data 



The MACHO Project ( lAlcock et al.l 119971 ) comprises eight years worth of observations 
of the Large and Small Magellanic Clouds and the B ulge of the Milky W ay. Our sample is 
drawn from the MACHO LMC Variable Star Catalog (lAlcock et al.ll2003l ). Sources from the 
full MACHO database of several million objects were selected for the Variable Star Catalog 
if the central 80 percent of points in the object's light curve failed to fit a constant magnitude 
in a x-squared test. This criterion resulted in 207,632 candidate variables in the LMC0 

MACHO data were taken simultaneously in two non-standard fil ters: Red and Blue. 
These can be transformed using the method of Alcock et al. Jl999bh to Cousins V and 
- R, and then used to find the Wesenheit reddening free magnitude, W = R — 4(V — R) 
( lAlcock et al.lll995l ). The construction of W allows very dim stars to enter our sample, 
so we have removed stars beyond MACHO's dim limit io Our LMC sample is defined by 
V — R > 0.5 and W < 15 (see Figured]) and is composed of 56,453 stars. These luminosity 
and color limits encompass 98 percent of the LPVs from Paper I. 

MA CHO employed a nonparametric phasing technique known as the SuperSmoother 
Method (jReimannl 1 1 9941 ) to phase all of the light curves in the Variable Star Catalog. This 
technique is robust against complex light curve morphology, but fails when presented with 
strongly multi-periodic behavior. Multi-periodicity is common among LPVs, and Super- 



1 The MACHO Project is available online at: http://wwwmacho.anu.edu.au/ 
2 Stars with instrumental magnitudes > —1.74. 
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Fig. 1. — The W vs. V — R color-magnitude diagram of the MACHO LMC Variable Sta r 
Catalog. W = R — 4(V — R) is the Wesenheit reddening free magn itude (lAlcock et al.lll995[) . 



The highlighted objects are those identified as LPVs in Paper I ( IFraser et al 



2005|), which 



used the SuperSmoother period and K s to determine variable type. Our sample is drawn 
from those stars with V — R > 0.5 and W < 15. This region surrounds the AGB and 
encompasses 98 percent of the LPVs from Paper I; the remaining 2 percent, which lie bluer 
than the present sample, are galactic foreground stars with long SuperSmoother periods. 
SuperSmoother failed to find a period for approximately half of the stars in this sample, 
while our technique succeeds 87 percent of the time. 
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Smoother failed to find a period for nearly half of our sample. 



We used the CLEANest algorithm of iFosterl (119951 . ll996a|jbT) to determin e the frequency 
ch aracteristics of the st ars in our sample. As implemented by lRorabeckl (119971 ). and described 



m 



Alcock et al.l (11999al ). CLE ANest uses the robust date-compensated discrete Fourier trans 



form (DCDFT) algorithm of iFerraz-Melld (119811 ). which finds accurate estimates of the 
amplitudes of the Fourier spectrum for data with uneven time sampling. The CLEANest al- 
gorithm iteratively finds the most significant peak in the power spectrum from the DCDFT, 
adds this frequency to those already known, determines the best fit by allowing all known 
frequencies to vary slightly, and subtracts the model light curve from the data. The algo- 
rithm exits when there is no longer any statistically significant powejfl in the frequency spec- 
trum. We veri fied our implementa tion by checking our results against th e 41 M ACHO Beat 
Cepheids from lAlcock et al.l (119951 ). and the test dataset from iRorabeck! (Il997l ). CLEANest 
successfully found every frequency down to our significance threshold. 

In the final analysis of our LMC sample, we searched a frequency space of 0.0003 day -1 
to 0.3 day -1 (corresponding to periods from 3.3 days to 3333 days) with a typical frequency 
resolution of 0.00003 day -1 . The average number of frequencies returned for a Blue light 
curve was 13. Although the MACHO light curves have the necessary time span to uncover 
frequencies as low as 0.0003 day -1 , the presence of power at these frequencies is likely caused 
by slow mean brightness changes. Our analysis found an excess number of frequencies 
between 0.0003 day -1 and 0.0006 day -1 (periods between 1666 days and 3333 days). We 
interpret this as being due to mean brightness changes over the timespan of the MACHO 
observations, and remove these frequencies from further analysis. 



We employ Two Micron All Sky Survey (2MASS — ICutri et al.l (120031 )) K s measurements 
as our luminosity indicator. The use of infrared luminosities splits Wood's Sequence B into 
our Sequences 2 and 3. We chose to take all 2MASS matches within two arc-seconds of the 
MACHO soured This is the distance at which there is a 50 percent chance of a false match. 
More than 95 percent of our sample has a match within this radius. 

Excluding light curves with fewer than 50 points (1 percent of our sample), and stars 
where the CLEANest analysis failed to converge to reasonable values (8 percent of our 
sample), the final number of stars with good Red and Blue CLEANest periods is 48,990 
(87 percent of our original sample). A comparison of the frequencies found in both the Red 



3 We use a limit o f 2, i n units of the DCDFT's power, based on the work of IRorabeck! <| 19971) and the 
suggestion in IFosterl (|l995l ) that power levels below 2 are "not even remotely significant." 



We used an updated astrometric solution for the LMC based on the UCAC system (jZacharias et al 



2000). 
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and Blue light curves shows that the error in our period estimates is within one-half of a 
percent up to periods of 55 days, after which it grows to approximately three times our 
typical frequency resolution (or 0.00009 day -1 ). 

The primary Fourier period of each star, in days, is the inverse of the first fundamental 
frequency found in the Blue light curve, Po = 1/vq. Plotting K s versus log 10 Po (insert 
of Figure [2]) immediately reveals the familiar period- luminosity sequenc es of the LPVs in 



the LMC. The sequences are named in the manner of iFraser et al.l (120051 ). from shortest to 
longest period: 4, 3, 2, 1, E, and D. The SuperSmoother and CLEANest analyses of our 
56,453 star sample is published in its entirety in the electronic edition of the Astronomical 
Journal. 



2.1. Artifact Removal From the Period-Luminosity Diagram 

A very strong vertical feature at periods of one year (log 10 P = 2.56) overlaps both 
Sequences 1 and D in the Fourier period-luminosity diagram, shown in the inse t of Figure [2] 



This fea ture is a result of the annual observing schedule of the MACHO Project (lAlcock et al. 



(jl999bl ). §6.2), and is not an intrinsic property of the star (a similar feature that corresponds 
to one month is faintly visible in Figure [2] at log 10 Po = 1-49). Since the strongest period 
of these stars is not due to the star itself, we associate these 11,215 stars with the One- 
Year Artifact rather than the sequences with which they overlap in period-luminosity space. 
Stars associated with the One- Year Artifact are, by their inclusion in the MACHO Variable 
Star Catalog, variable objects, but they are not necessarily periodic. They certainly have 
no higher amplitude periodicity in their frequency spectra than the weak signal due to the 
annual schedule of earthbound telescopes. 

In Paper I we simply masked all the stars in this region. In this work we identify just 
the stars associated with the One- Year Artifact by exploiting differences in the properties 
of the stars belonging to the artifact with the stars nearby in the Fourier period-luminosity 
diagram. We identify stars as associated with the One- Year Artifact if they lie in the region 
2.53 > log 10 Po > 2.60, and according to the following rules: 



• Sequence 1 is composed of the reddest stars in our sample. We identify stars as part 
of the One- Year Artifact if their luminosity places them in or above Sequence 1 and 
they are bluer than J — K s = 1.5. 

• At luminosities dimmer than Sequence 1 we find that the One- Year Artifact stars have 
poorly correlated Red and Blue light curves. We identify stars as part of the One- 
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log Period 



Fig. 2. — Fourier period-luminosity diagram of the LPVs (Long-Period Variables) in the 
Large Magellanic Cloud with and without the stars identified with the One- Year Artifact. 
Contour levels indicate the density of stars per 0.05 in log 10 Po and 0.1 in K s . The Fourier 
periods shown for stars in Sequence E are half of the orbital periods of these binary systems. 
When Sequence E is plotted at the orbital period (+0.30 in log space) it smoothly joins 
with the bottom of Sequence D. The inset in this figure shows the Fourier period-luminosity 
diagram with the One- Year Artifact stars included; the contour levels in the inset start at 5 
stars per 0.05 in log 10 P and 0.1 in K s , and continue at the same contour levels as the larger 
diagram. 
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Year Artifact if their Red and Blue periods differ by more than 20 percent, or if the 
amplitudes corresponding to those periods differ by more than 90 percent. 

• In the luminosity range between Sequences 1 and D we find that we require the ad- 
ditional parameter of the \ 2 statistic for a sine wave corresponding to Pq, which is in 
the range 2.53 > log 10 Po > 2.60. We select stars with either Blue amplitudes of less 
than 0.05 magnitudes mean-to-peak, or those with x 2 < 7 x 10 4 . 

Unfortunately the stars in the One- Year Artifact that are dimmer than Sequence D are 
difficult to differentiate from the background, although certainly the great majority of stars 
in this region should be identified with the One- Year Artifact. We chose to simply identify 
stars as members of the One- Year Artifact if they lie below Sequence D and are in the normal 
period range of 2.53 > log 10 P > 2.60. 

The One- Year Artifact consists of 11,215 stars with very weak or nonexistent periodicity, 
or 24 percent of all the stars with good Red and Blue CLEANest results and a 2MASS K s 
magnitude. This greatly outnumbers other stars in this narrow period range. Although our 
process likely does not isolate every star that is associated with the One- Year Artifact, it 
does serve to uncover what is present below this distracting feature. The Fourier period- 
luminosity diagram of the LMC with the One- Year Artifact removed is presented in the main 
panel of Figure [2j 



2.2. Finding the Average Infrared Light Curves 

Although the amplitude of LPVs in the infrared is much lower than in the visibl^l, 
there is still intrinsic scatter in the 2MASS observations due to the measurement of each 
star at a random point in its light curve. C ollections of stars wit h similar light curves can be 



corrected for this effect in the manner that iNikolaev et al.l (120041 ) used for Cepheid variables. 



The Fourier period-luminosity sequence in some band is fit with a relation that includes a 
correction term which is a function of <ft, the phase of the 2MASS observation with respect 
to the MACHO light curve. So for each star, i: 

mi = a\og l0 Pi + [3 + {}{&) (1) 



5 In classical pulsators this is due to the black-body behavior of the atmosphere's continuum emission. 
The large amplitudes of Mira light curves cannot be modeled by this behavior alone; it is also necessary to 
consider very strong effe cts from the formation of molecules in the photosphere of the star. See, for example, 



Reid k Goldstonl (|2002h 
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Fig. 3. — Example optical and infrared light curves for the star with MACHO designation 
55.3126.13. The infrared light curve is taken from the 2MASS calibration tile 90400, and 
both are phased to the CLEANest period. The infrared light curve shows a phase lag of 
approximately 10 percent as compared to the optical (the vertical dashed line indicates the 
approximate maximum of the infrared light curve). This star has an optical period of 64.03 
days and lies in Sequence 3. It's optical light curve is shown with the CLEANest model of 
the primary period; all other frequencies have been subtracted from the model. The average 
infrared light curve for Sequence 3 was too low in amplitude to measure using our technique, 
so this star must have unusually strong variation. 
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The correction term, fl(4>i), is the average infrared light curve for these stars. We chose 
a form for the correction function based on 2MASS light curves available for 46 stars from 
our sample. Light curves exist for stars in the 2MASS "calibration tiles", fields which were 
observed multiple times each night to provide photometric calibration. Several tiles lie in 
the vicinity of the Magellanic CloudsJ^I but only tile 90400 provided light curves of sufficient 
length to investigate the behavior typical of long-period variables. An example 2MASS light 
curve phased to the period of the corresponding MACHO Blue light curve is shown in Figure 
[3J For the 46 stars from our sample that match 2MASS sources from this tile, we found that 
a second-order Fourier series was an adequate model of the light curves in J, H, and K s . 



ft(0) = ^ Aj cos(2ttj0) + Bj sin(27rj0) (2) 

The phase of the 2MASS observation, <pi, for each star i, is the fractional part of the 
difference between the time of maximum light, T max j, for the primary Fourier period and 
T 2mass i, the time of the 2MASS observation, in units of the period. 



modf 1 (3) 



The model of the infrared light curve, as based on the 46 2MASS light curves, is fit si- 
multaneously with each period-luminosity relationship in our f ull sample . The scatter about 
each of these relationships is not completely accounted for, e.g. Ilta et al.l (120041 ) found differ- 
ent relations for RGB and AGB stars. For this reason we hav e fit each sequence separatel y 
above and below the tip of the Red Giant Branch (K s = 12.3 (iNikolaev fc Weinberg l2000l )) 
for Sequences 2, 3, and 4. Additionally, the widely varying properties of the stars that 
are being combined tends to reduce the amplitude of the average light curve. We take the 
limit of detectability of the infrared light curves to be a mean-to-peak amplitude of 0.02 
magnitudes — the typical magnitude error of our 2MASS observations — and further require 
that the scatter about the period-luminosity relationship narrows with the addition of the 
infrared light curve. At this level we detect infrared light curves in J, H, and K s for Se- 



6 These tiles, 90298, 90299 in the SMC and 90400, 90401, and 90402 in the LMC, were defined to support 
the 2MASS deep observation campaign of the Magellanic Clouds. Of these five tiles, and the calibration tiles 
from the rest of the survey, only 90400 and 90401 overlap with MACHO observations. We only used data 
from tile 90400 because it was visited more often (377 times, versus 156) over more nights (approximately 
40, versus approximately 13) and over a longer time span (90 nights, versus approximately 30) than tile 
90401. 
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quence 2 above the tip of the RGB, and in Sequences 1, E and D. The fitting and correction 
was performed individually for two groups in Sequence D: those with Blue P peak-to-peak 
amplitudes less than or equal to 0.2, and those with amplitudes greater than 0.2. This 



separat ion roughly corresponds to the division of LMC stars in OGLE by ISoszynski et al. 



( 12004al ) into OSARGs and SRV/Miras. In Sequence 1 we found it necessary to create three 
amplitude bins (with boundaries at 0.4 and 1) as well as separating the oxygen and carbon 
stars (using the color cut J — K s = 1.4). This created a total of six groups in Sequence 1 
that were each fit individually. The results of these fits are tabulated in Table [TJ and the 
infrared light curves are discussed in §3.11 



In lNikolaev et al.l (120041 ) the remaining scatter around the period- luminosity relationship 
for LMC Cepheids was used to derive individual distances and extinctions for each star. In 
that work, performing the original fit again after fitting for the distance and extinction 
together yielded a closer fit in all bands. We have tried this technique in J, H, and K s using 
the scatter around the Fourier period-luminosity relationship in W as a test of its success 
(W cannot be used to constrain extinction since it is constructed to remove the effects of 
reddening). Unfortunately the dispersion in W failed to improve after this step, leading 
us to believe that the extinction in these stars requires a more careful model than simply 
relying upon the color excess. In a separate test, we found that the addition of a color term 
also failed to improve the dispersion in the unrelated Fourier period-luminosity relations. 
LPVs must have some intrinsic luminosity variation that is not well modeled simply by 
using pulsation with color or individual distances and extinctions. 

Our infrared light curves have very low amplitudes due to the varying properties of their 
constituent stars. Because the correction is so small, any individual star's departure from 
the average light curve is likely to dominate the effect of the random phase of the 2MASS 
observation, so we have not applied these corrections to the 2MASS magnitudes used in this 
work. While the improvement of the Fourier period- luminosity fits is modest, the infrared 
light curves do prove useful in providing physical insight into the variability mechanisms as 
described in £13.11 below. 



3. Results 

The Fourier period-luminosity diagram of the LMC with the One- Year Artifact removed 



i s prese nted in Figure [2J The sparsely populated fifth sequence found by ISoszynski et al. 



( I2004al ). which lies toward shorter periods of Sequence 4, is not seen in this diagram. This 
sequence is only seen in the secondary periods (-Pi), and we have plotted just the primary 
Fourier period (Pq) for each star. All of the stars that compose the fifth sequence have 



13 



stronger periods elsewhere on the Fourier period-luminosity diagram, primarily on Sequence 
4 or the One- Year Artifact. 

Stars in this diagram were grouped into sequences based on the contour traced at a 
density of five stars per 0.05 in log 10 Po an d 0.1 in K s (as shown in t he inset of Fig ure El) . 



The stars in the high-luminosity tip of Sequence 1 are carbon stars (jGroenewegenl 1200 



some of which are heavily self-extincted and fall into the gap between Sequences 1 and D . 



Such stars are easy to recognize due to their heavy reddening (jNikolaev fc Weinberg! 120001 ) . 
thus stars in the gap between Sequences 1 and D with J — K s > 1.4 are assigned to Sequence 
1. 

Se quence E is well known to merge with Sequ ence D when plotted at its true orbital 



period (jSoszynski et al.ll2004bl ; iDerekas et al.ll2006l ). but we have chosen to plot Sequence 



E at its Fourier period like the other stars on this diagram. This allows us to differentiate 
these stars from those at the bottom of Sequence D, but it presents us with the problem of 
separating the two sequences. It's clear that Sequence E in Figure [2]is poorly populated, and 
at the top of the sequence (K s < 13) it's difficult to tell where Sequence E exists among the 
background of stars between Sequences 1 and D. We chose limits between these sequences 
such that Sequence E would remain sparsely populated at its high-luminosity end. 

Compared to Paper I, where half of that sample had no period assigned by the Super- 
Smoother analysis, or a period assigned to one year or multiples of one day, we have added 
approximately 20,000 stars to the period-luminosity diagram. The newly analyzed stars lie 
mostly on the sequences with the lowe st amplitudes (D, 4, and to a lesser extent, 3) and are 



below the tip of the RGB (K s = 12.3, iNikolaev fc Weinberg! (120001 ) ) . An example is shown 



in Figure @J top panel. Sequence E stars were well represented in Paper I even though their 
light curves also have low amplitudes. Half as many stars are identified with the One- Year 
Artifact due to the careful identification of these stars, as opposed to masking all stars with 
in this region. Luminosity functions for many of the sequences are shown in Figure [51 

Sequence D as presented in Paper I was underrepresented, accounting for only 9 percent 
of LPVs. It now represents 31 percent of the stars with good Red and Blue CLEANest 
results and a 2MASS K s magnitude (Table [2]). The right panel of Figure El shows the close 
match between the luminosity functions of Sequence D and our sample, suggesting that stars 
on Sequence D are drawn from the entire population of long-period variables in the LMC 

The amplitudes of the light curves that are quoted in this work are the peak-to-peak 
amplitudes found by CLEANest for the Pq term. CLEANest amplitudes underestimate the 
actual light variation since some of the power associated with this period is contained in 
harmonics and mixing terms. The amount by which the amplitudes are underestimated is 
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Fig. 4. — Example light curves and CLEANest models of LPVs from our sample; each panel 
is labeled with the star's MACHO designation. The top panel shows a star from Sequence 
4 with close Fourier periods and a sufficiently complicated light curve that a period close to 
1 day was falsely assigned in Paper I. The second panel shows another star from Sequence 4 
that exhibits amplitude modulation due to its beating periods of 52.18 and 53.40 days. The 
bottom panel shows the entire MACHO light curve of a Mira variable from Sequence 1. This 
star's periodicity changes character over the course of the survey, and it has a CLEANest 
spectrum with beating periods of 168.55 and 169.58 days. A representative sine curve at the 
primary CLEANest period is shown. 
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Fig. 5. — Luminosity functions of selected groups of stars from our sample. The panel 
on the far left shows the luminosity function of the One- Year Artifact. The second panel 
compares Sequence 4 (black outlined bars) with Sequence 3 (gray bars with no outlines). 
The third panel similarly compares Sequences 2 (black outlined bars) and Sequence 1 (gray 
bars with no outlines). Finally, the panel on the far right compares the luminosity functions 
of Sequence D (black outlined bars) with our sample (gray bars with no outlines, scaled by 
1/3). 
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found by comparing the CLEANest amplitudes to amplitudes listed in the MACHO Vari- 
able Star Catalog for stars with SuperSmoother periods. SuperSmoother amplitudes were 
calculated by finding the difference between the mean magnitudes of the closest points to 
the maximum and minimum of the SuperSmoother phased light curve. On average, the ratio 
of the SuperSmoother amplitude to the CLEANest amplitude is 1.5 for stars in Sequence E 
and stars in the long-period edge of Sequence 1, while the remaining stars on the Fourier 
period-luminosity diagram tend to show ratios of 1.7. We have not applied a correction 
factor to our CLEANest amplitudes, and use them only for relative comparisons. 

The relationship between the pulsation amplitude and log 10 Po is markedly different for 
the different sequences. Figure |6] shows that for stars in Sequences 1-4 there is a correlation 
between increasing period, increasing amplitude, and increasing range in amplitudes. As 
compared to stars on Sequences 3 and 4, stars on Sequences 1 and 2 pulsate with much higher 
amplitudes. Sequence E, as expected for binary stars, does not show a strong dependence 
of amplitude on period. This sequence is a clear continuation of the lower amplitude group 
of Sequence D stars, those related to the OS ARGs. The modera te amplitude-luminosity 



correlation for stars in Sequence D reported by iDerekas et al.l (120061 ) is not seen for the bulk 



of the small- amplitude population of Sequence D, using log 10 P as a proxy for luminosity. 



3.1. The Average Infrared Light Curves of LPVs 

The comparison of the average infrared light curves with the corresponding optical light 
curve's properties is a useful constraint on variability mechanisms in these stars. In many 
types of pulsating stars — including Miras — pulsation has a stronger effect in the optical than 
the infrared, while variation due to binary systems show light curves of similar amplitudes 
in all wavebands. Although the average amplitudes of infrared light curves are smaller than 
the specific star's amplitudes that they are fit to, comparisons can be made in a relative 
sense. While Table [1] presents information for all of the infrared light curves, a summary 
of amplitude ratios and phase lags for the average Blue and K s light curves is presented in 
Table [2] for each of the sequences. 

We find a clear division in the B\ue/K s amplitude ratios between Sequences 1 and 2, 
and Sequence E. As expected for binary systems, Sequence E has a smaller amplitude ratio 
and it has similar amplitudes among the J, H, and K s light curves. The pulsating stars 
in Sequences 1 and 2 show higher ratios and decreasing amplitude with redder wavebands. 
Interestingly, only the low-amplitude group of Sequence D stars is similar to Sequence E. 
Although the higher amplitude group in Sequence D has similar amplitudes among the 
2MASS bands, the average Blue amplitude is much higher, leading to an amplitude ratio of 
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log Period 



Fig. 6. — Peak-to-peak MACHO Blue amplitude versus log 10 P f° r stars on Sequences 1-4, 
and E and D. Stars are color-coded according to the sequence on which their primary Fourier 
period lies, as shown in the inset. The Fourier periods shown for stars in Sequence E are 
half of the orbital periods of these binary systems, a difference of +0.30 in log space. 
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eight, which is most similar to Sequence 1. 

Pulsation modes in LPVs c an also be constrain ed using the phase lag between the optical 



and the infrared light curves. ISmith et all (120061 ) searched for such phase lags using data 



from the DIRBE instrument on the COBE satellite. In their sample of 21 stars, all of the 
Miras — including one carbon star — showed phase lags, while four of the five SR variables did 
not. They compared the stars in their sample to time- resolved dynamical models of oxygen- 
rich stars from the literature, and found that phase lags are predicted for fundamental-mode 
oxygen rich stars, but not for stars pulsating in the first overtone mode. The carbon star 
models available at the time did not consistently predict a phase lag. Our average infrared 
light curves show phase lags of 10-20 percent for both Miras and SR variables. The example 
light curve from the 2MASS calibration tiles, shown in Figure [31 demonstrates this phase 
lag. The average infrared light curves for the 0.4 < amplitude < 1 bins in Sequence 1 are 
the exceptions to this rule. Neither the oxygen stars nor the carbon stars in this range show 
the typical phase lag commonly observed among other LPVs. However, the most notable 
exception is Sequence D, both its low-amplitude and high- amplitude infrared light curves 
lead the optical light curves, unlike any of the other sequences. 



3.2. Multi-Periodic Stars 



Stars on Sequence D are well-known to be multi-periodic pulsators; IWood et al.l (119991 ) 
found that Sequence D stars exhibited a shorter period that fell on his Sequence B (which 
is composed of our Sequences 2 and 3). Figure [7] shows the secondary Fourier period (Pi) of 
all of the Sequence D stars overlaid on the normal Fourier period-luminosity d iagram. We 



see th at many of these periods do in fact lie on Sequences 2 and 3 as found by IWood et al. 



(119991 ). although we also see that Sequence 3 is favored. There are also a substantial number 
of secondary periods of Sequence D stars on Sequence 4 below the tip of the Red Giant 
Branch, and some between Sequences 1 and 2. However, many of these secondary periods 
still fall on Sequence D, or have periods half as long as in Sequence D. These stars with 
Pq/P\ ~ 1 and Pq/P\ ~ 2 are discussed below in the context of the Petersen diagram. 

The Petersen diagram is the plot of the ratio of the two strongest periods versus the 
longer period in a multi-periodic star, and is presented in Figure Multiple periods are 
found for the great majority of stars in our sample, so stars from all of the Fourier period- 
luminosity sequences are represented here. The One- Year Artifact is visible in this plot since 
we only removed stars whose primary Fourier periods have low or nonexistent periodicity. 
Stars with their secondary Fourier periods on the One- Year Artifact lie either in a vertical 
stripe at P = log 10 (365 days) (when their secondary period is longer than the primary, since 
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Fig. 7. — Fourier Period-luminosity diagram showing secondary Fourier periods (Pi) for stars 
with their primary Fourier period (P ) on Sequence D. Light gray points show the entire 
sample of P Fourier period-luminosity relations. 
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log Long Period 

Fig. 8. — Petersen diagram of the LMC: for the primary and secondary Fourier periods of 
a star, the ratio of the longer to the shorter versus the log of the longer period. Stars are 
color-coded according to the sequence on which their primary Fourier period lies, as shown 
in the inset. The structure in this diagram is discussed in §3.21 
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it is the longest period that defines pla cement on t h e hor izontal axis) or along the curve 



Pol Pi — -Po/(365 days). As observed in I Wood et al.l (119991 ). there is a wide locus of points 



centered around period ratios of 10; these represent stars with their primary or secondary 
Fourier periods lying on Sequence D. 

There is a great deal of structure lying at period ratios below three, which we show in 
a variant of the Petersen diagram: Figure [9] plots the Pq/P\ ratio for our stars with respect 
to the log 10 P Fourier period. Note that P is the strongest period, but not necessarily 
the longest period, as shown by the existence of many stars with period ratios of less than 
one. In this plot, stars with Pi on the One- Year Artifact lie only along the curve Pq/P\ = 
Po/(365 days). Aside from stars with periods on the One- Year Artifact, there are at least 
10 groups of stars visible with period ratios between zero and two. 

The groups of stars with the smallest period ratios, and therefore the longest secondary 
periods (Pi) relative to their primary period (Po), have their primary periods lying on one 
of the numbered sequences and their secondary periods lying on Sequence D. 

Many stars show period ratios of approximately one. These closely spaced frequencies 
are often taken to indicate non-radial pulsation among the least luminous stars from all of 



the n umbered sequences; they ha ve been seen before in b oth Galactic field stars (IKiss et al. 



20001 ) and the Magellanic Clouds (jSoszynski et al.ll2004al ) . The solid lines indicate the region 



where the period ratio is exactly one within our estimated errors; arcs that are visible at the 
long-period extreme of this region are due to the limits of our frequency resolution. 

These frequencies can be seen to produce beats, as seen in the light curve of a Sequence 4 
star in Figure IU middle panel. It is also possible that some of these closel y spaced frequencies 



may be due to a single, slowing varying period. iTempleton et al.l ( 120051 ) found that approx- 
imately one-tenth of their sample of 547 AAVSO Mir a light curves showed "meandering" 
periods that did not simply increase or decrease in length. An example Mira variable from 
Sequence 1 that has closely spaced frequencies in its CLEANest spectrum, and a changing 
period, is shown in Figure HI bottom panel. It is overlaid with a sine curve corresponding 
to the primary Fourier period that CLEANest found: 168.55 days. If the light curve of this 
star is split in half at the 1760 day mark, we find that the Fourier period of the first half is 
166.62 days while the second half has a Fourier period of 163.86 days. 

Apart from the period ratios above, stars in Sequences 1 and 2 and stars in Sequences 3 
and 4 tend to show secondary periods that correspond to the other Fourier period-luminosity 
sequence in each group. The Sequence 2 stars at a period ratio of 0.6 have their secondary 
periods on Sequence 1, while the Sequence 1 stars at ratios between 1.6 and 2.2 have their 
secondary periods on Sequence 2. Likewise, stars on Sequences 3 and 4 show groups at period 
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log Primary Period 

Fig. 9. — Pq/P\ versus log 10 P showing just the period ratios below 3. Stars are color-coded 
according to the sequence on which their primary Fourier period (Pq) lies, as shown in the 
inset. Note that Pq is the strongest period, but not necessarily the longest period, as shown 
by the existence of many stars with period ratios between and 1. Stars with secondary 
Fourier periods (Pi) on the One- Year Artifact lie along the curve Pq/P\ = Po/(365 days). 
The structure in this diagram is discussed in §3.21 
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ratios of 1.4 and 0.7. We do see stars whose periods cross over between these two groups, 
for example, the Sequence 4 stars at a period ratio of 0.5, which indicate a secondary period 
on Sequence 2, and a few Sequence 3 stars at a period ratio of 0.65, indicating a secondary 
period on Sequence 1. There are even fewer stars that cross from Sequences 1 or 2 to the 
shorter period sequences, but some Sequence 2 stars are visible at period ratios around two, 
which represent secondary periods on Sequence 4, and a scattering of Sequence 1 stars lie at 
period ratios greater than three, which represent secondary periods on Sequence 3. The fact 
that the period ratios of stars on each pair of Sequences, 3 and 4, and 1 and 2, are reciprocals 
means that on a Petersen diagram these groups would substantially overlap, indicating that 
similar physical mechanisms are occurring in each case. 

Visible at a period ratio of exactly 0.5, 1.5 and 2 are groups of stars from Sequence E. 
The CLEANest periods found for these contact binaries represent only half of the orbital 
period of these systems; therefore, these period ratios represent the true orbital period, and 
the third and fourth harmonics, respectively. These stars often have minima of alternating 
depth which is characteristic of contact binaries (Figure [TO]) . 

The Sequence D stars with Pq/P\ ~ 2 have period ratios of exactly two within our 
estimated errors, indicating the prevalence of the second harmonic in these light curves. The 
cause for the closely spaced frequencies (Pq/P\ ~ 1) in Sequence D stars is still unknown at 
present. 



Finally, the sparsely populated fifth sequence (jSoszynski et al.l l2004al ) is visible in the 
Sequence 4 stars with ratios of 1.3. 



3.3. Period Changes 



P eriod changes are well known for Miras (ITempleton et al.ll2005l ) and semi-regular vari- 
ables (IKiss et al.l |2000| ). iGroenewegenl (120041 ) observed changes from historical periods in 
OGLE LMC LPVs that indicated stars had moved between Sequences 1, 2, and E0. It's 
likely that some of the stars in our sample underwent a period change during the eight years 
of the MACHO Project. To explore this possibility, we have split each of our light curves into 
two halves of equal length and run our analysis on each half independently. The split light 
curves are almost four years long, which lLah et al.l (120051 ) found was a sufficient timespan to 
resolve the familiar Fourier period-luminosity sequences. This is the case for our split light 



7 Actually between Sequences B, C, and D. We take their Sequence B as S equence 2, but not S equence 3, 
because these stars lie on the long period side of Sequence B (Figure 8 from IGroenewegenl (|2004|)). 
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Fig. 10. — Example light curve of a star on Sequence E with the characteristic alternating 
minima of a binary system. Note that this light curve is phased at twice Pq period. 
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curves as well, but instead of finding changing periods, this technique appears to find stars 
with multiple periods of almost equal power. 

Fifty percent of this sample shows changes in log 10 period of greater than 0.1 dex, the 
approximate width of the sequences. We have tabulated these movements in Table [2] for all 
cases where more than 10 percent of the stars in one sequence move to another. Many move 
to Sequence D, and we also see stars switching sequences between the pairs 3 and 4, and 1 
and 2. These movements between the pairs of numbered sequences only happen for stars in 
the inside edge of each pair of sequences, as shown in Figure [HQ 

Considering the multi-periodic nature of these stars, a nd the much smaller observed 



rate of period changes in Miras from iTempleton et al.l (120051 ). we do not believe these data 



indicate a flurry of period changes. Instead, we observe that the changes we see correspond 
very well to the multiple periods explored in §3.21 Examination of some of the light curves 
of these stars shows that both periods are typically present in the frequency spectrum of 
both halves of the light curve, but the relative amplitudes of these periods change. For stars 
moving between a short period sequence (such as 3 or 4) and Sequence D we see that the 
power in the high frequency component becomes distributed among multiple closely-spaced 
frequencies, which no longer have greater power (individually) than the stable low-frequency 
component. It's unclear if this is an effect of the different time-sampling in each half of the 
light curve, or if it is evidence of some change in the star itself. 

Although this analysis may not find period changes directly, th e equivalence of these 



periods can be understoo d as a result of longer term period changes. IWhitelockl (119861 ) and 



Lattanzio fc Wood! (120031 ) argued that the period-luminosity sequences could be understood 
as the result of different pulsation modes that are excited in turn as the star evolves. The 
stars that are close to the point where the dominant modes switch (and therefore the sequence 
switches as well) would be expected to have multiple periods with similar amplitudes 



4. Discussion 



The analysis of the frequency spectra of variable stars has been used with great success 
for characterizing light curve morphology, identifying binary s tars, and constrai ning observed 



pulsa tion modes in studies of Cepheids (ISimon fc Led 1 1981), Be at Cepheids ( Alcock et al. 



1998|), RR Lyrae (lAlcock et al. 



1995), Type II Cepheids and RV Tauri stars ( Alcock et al.1 

200oh . and LPVs (see JQ). Here we discuss how these techniques contribute to the under- 



standing of the origin of Sequence D and the relationship between the different LPV stages 
and stellar evolution. 
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log Period 



Fig. 11. — Fourier period-luminosity diagram showing period changes of less than 0.5 dex, 
color-coded by the log 10 change in period between the first half of the light curve and the 
second half. Stars are plotted at the period found in the first half of the light curve. The 
majority of the period changes of stars in Sequences 4, 3, and 2 are to Sequence D (see 
Table EJ); these larger period changes (which are typically > 0.5 dex) are not plotted so 
that the underlying pattern of smaller period changes is visible. These smaller changes show 
switching between Sequences 3 & 4, and Sequences 1 & 2, but only for stars on the inside 
edges of these pairs. 
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4.1. Sequence D 



One-third of the stars in our sample exhibit the "long secondary period" phenomenon. 
Although the exact mechanism for Sequence D is still unknown, there are many reasons to 
think that it is correlated with binary systems. In this paper, the evidence includes the 
similarity between the luminosity functions of Sequence D and our entire sample (Figure 
EJ right panel), and the similar amplitudes of Sequence D's average infrared light curves 
across the 2MASS bands. Other evidence includes the smooth connection between Seq uence 
D and Sequence E — when E is pl otted at its orb ital period (jSoszynski et al.l l2004bl ) . the 
presence of ellipsoid al light curves (lSoszyriskill2007l ) . and radial velocity studies of these stars 
( lAdams et al.ll2006l ). However, there are also several significant features of Sequence D that 
are unique, or show differences from Sequence E. Sequence D does not suffer from period 
halving like Sequence E does, and Sequence E stars often have the third and fourth harmonics 
present in their frequency sp ectra, while Sequence D stars show the second harmonic instead 
( §3.2p . iDerekas et al.l (120061 ) found that in an amplitude-luminosity plot, stars in Sequence 
D follow a different pattern than in Sequence E. We see that Sequence E appears to be 
a continuation of only the small- amplitude (< 0.2 Blue peak-to -peak) group of stars i n 
Sequence D, those stars that roughly correspond to the OSARGs of ISoszynski et al.l (j2004al ). 
This is consistent with the comparison of the average light curve amplitudes in Blue and 
K s , where only the lower amplitude Sequence D stars were a good match to Sequence E. 
Finally, Sequence D's infrared light curves lead the optical light curves by 10-15 percent, 
which is a feature unique to this sequence. These facts do not necessarily preclude a binary 
star mechanism for Sequence D, but they are useful constraints for proposed mechanisms. 
We note, as an additional constraint, that stars associated with Sequence 1 are far less likely 
to have a period lying on Sequence D, and that Sequence D is not observed in the LMC 
below K s « 13.7. 



The model pr oposed by ISoszyhskil (120071 ). based on the original model proposed by 
Wood et al.l (119991 ). is that of a binary system in which the mass lost from the red giant is 
concentrated near the companion, and regularly obscures the red giant. The wide range of 
observed light curve amplitudes for Seq uence D stars, from 0.1 up to five magnitudes in the 
MACHO Blue filter (TFraser et alJliooj ). can be readily explained by the projection effects 
of different inclination angles in a binary system. 

If Sequences E and D are truly composed of binary systems, then the population of 
binary stars in our sample includes stars with either their primary or secondary Fourier 
period lying within the boundaries of these sequences (e.g. the top panel of Figure 0] shows a 
star whose secondary Fourier period lies on Sequence D). After removing periods identified 
with the One- Year Artifact, we find that 48 percent of the stars in our sample have variability 



-28- 



associated with Sequences E or D (see Table [2]). For comparison, iDuquennoy fc Mayor 
(199 1) found a bina ry fraction of approximately 40 percent for nearby solar-type stars, and 



Reid &: Gizisl (119971 ) found a fraction of a pproximately 35 percent for low-mass stars, a trend 



in mass which is discussed in lLadal (120061 ). It is reasonable to assume that we cannot see pole- 
on binaries, and there is no reason to think that all binaries have periods shorter than four 
years, both of which imply that 48 percent underestimates the total percentage of binaries 
seen in the LMC This is a serious problem for any explanation of Sequence D that relies 
solely on binary systems. However, it is very likely that a subset of these stars do show 
variability due to binarity, perhaps stars in one of the populations that can be separated by 
color or amplitude. 



4.2. Comparison to Evolutionary Models 



We can begin to characterize the e volution of LPVs by comparison in color-magnitude 
space to models. iMarigo et al.l (120031 ) produced a population synthesis model of the K s 
vs. J — K s color-magnitude diag ram (CMP) o f the L MC using the RGB and Early AGB 
(E-AGB) evolutionary models of iGirardi et al.l (120001 ). and a preliminary version of their 
thermally pulsing AGB star models (published later in IMarigo &: Girardil (120071 )). Their 
Figure 12 illustrates the luminosities and colors corresponding to major phases in an LMC 
star's evolution from RGB, through the Early AGB, and finally to the thermally pulsing 
AGB (including a transformation to carbon-dominated atmospheres for some stars). For 
comparison, our Figure [12] shows the average J — K s color binned in magnitude for each of 
the Sequences 1-4, E, D, as well as the One- Year Artifact and the background population 
of stars. The color of the "background" stars match Galactic disk t urn-off stars a nd LM C 
intermediate-mass stars on the Early AGB in the synthetic CMD of IMarigo et al.l (120031 ). 



As shown in Figure 12, the J — K s color is more effective at distinguishing evolution- 
ary stages in the AGB than in the RGB. However, with the additional information from 
the luminosity functions, we can also estimate the importance of RGB stars to each se- 
quence. A distinct peak in the luminosity function at the tip of the RGB (K s = 12.3, 
Nikolaev fc Weinberg (120001 )) is widely taken to indicate that the majority of the stars dim- 



Ita et al. 


2004; 


Fraser et al. 


2005) 



The stars in our sample that show very weak or nonexistent periodicity (the 24 percent 
identified with the One- Year Artifact) are predominantly RGB stars. The left two panels 
of Figure [5] show the luminosity functions of stars in the One- Year Artifact with Sequences 
3 and 4. Below the tip of the RGB the population of One- Year Artifact stars dominates, 
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Fig. 12. — Average J — K s color for each sequence, as well as stars in the One- Year Artifact 
and stars in the "background" of the Fourier period-luminosity diagram. J — K s colors are 
calculated for 0.5 magnitudes bins and shown for bins with more than 40 stars. 
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suggesting that very weak or non-periodic variability is common among RGB stars. Above 
the tip of the RGB there is a much closer correspondence between the One- Year Artifact and 
Sequences 3 and 4. Thus it appears that stars at the dimmest luminosities in our sample 
vary aperiodically while on the RGB, but most begin to show periodic behavior when they 
brighten to K s pa 13. 

After passing off of the tip of the RGB, stars may pulsate with shorter periods and lower 
luminosity as RR Lyrae on the horizontal branch. They become LPVs agai n as they ascend 
the Ea rly AGB (i.e. prior to the first thermal pulse or helium shell flash). ISoszynski et al. 



( I2004al ) used the slight offset in period between RGB and AGB stars to show that AGB 
stars pulsate alongside RGB stars below the tip of the RGB. 

Evolution proceeds to brighter luminosi ties until the onset o f thermal pulses, which be- 
gin at K s pa 12 on the synthetic CMD fr om iMarigo et al. rl2003h. Stars primarily populate 
Sequences 1 and 2 above this luminosity. IWagenhuber fc Tuchmanl (119961 ) predict that ther- 
mal pulses will create large modulations in luminosity and the pul sation period (and mode) 
in AGB stars with timescales of thousands of years. The models of IMarigo fc Girardil (120071 ) 
substantially agree with these predictions, and show large changes in pulsation period due 
to mod e switching as a direct result of a thermal pulse. Period changes in LPVs are well 
known (jTempleton et al.ll2005l ) but only a small percentage of stars at any one time should 
be undergoing a thermal pulse due to the short timescales of thermal pulses relative to the 
long inter-pulse period. The observed period changes of LPVs are not well explained by the 
effects of thermal pulses alone. 

The middle two panels of Figure [5] compare the luminosity functions of the numbered 
sequences. The relative importance of the two giant branches shifts from the RGB to the 
AGB as we move to Sequences 1 and 2. Additionally, the peak number of stars above the 
tip of the RGB in each sequence is found at highe r luminosity from Seque nce 4 to Sequence 
1. The OSARG versus SRV/Mira distinction of ISoszynski et al.l ( )2004al ). by virtue of its 
definition, roughly corresponds to a division between two pairs of sequences: Sequences 3 
and 4, and Sequences 1 and 2. This division is als o clearly seen i n the observed period ratios 
(Figure E]). Considering the synthetic CMD from IMarigo et al.l (120031 ) . OSARGS are closely 
related to RGB and E-AGB stars, while the SRV/Mira stars are more closely related to 
thermally pulsing AGB stars. 

At K s pa 11 the synthetic CMD predicts the formation of carbon stars, and the typical 
J — K s colors of each sequence diverge (Figure IT21). Only Sequences 1, 2 , and D redden 
to the expected J — K s color of the carbon star tail of IMarigo et al.l (120031 ). Their models 
also show that these stars do not evolve in brightness after this stage, so the observed range 
of carbon star luminosities in our sample may be interpreted as a range of stellar masses. 



31 



Marigo &: Girardil (120071 ) show that only stars between 1.2 and 2.5 M & undergo a dredge-up 
that can bring the results of nuclear burning to the atmosphere without quickly destroying 
it through hot botto m burning. Since M iras exist on Sequence 1 at luminosities both above 
and below K s = 11 (IFraser et al.ll2005l ). not all large a mplitude pulsators ar e carbon stars. 
Also, not all carbon stars have such red J — K s colors: iGroenewegenl (120041 ) found carbon 
stars on the shorter period sequences — the popular color-cut of J — K s > 1.4 appears more 
effective at segregating M stars, which are rarely this red. 

Using J — K s > 1.4 to select areas that include only carbon stars, we see that many 
occup y the highest luminosities of both Sequences 1 and 2, as also seen in lLebzelter fc Wood 
(120071 ). App roximately 40 perce nt of the carbon stars lie on Sequence 2, similar to the 
prediction of iMarigo et al.l (120031 ) . who fit the observed K s vs. J — K s CMD by assuming 
a 50 percent mix of fundamental and first-overtone pulsation among the carbon stars in 
their model. Figure [12] shows that the stars on Sequence 1 evolve to redder J — K s colors 
than on S equence 2, presumably due to increased mass loss driven by fundamental-mode 
pulsation (IMarigo fc Girardil 120071 ). Some stars on the long-period extreme of Sequence 1 
are underluminous for their color, which may be self-extinction due to dusty outflows. 

Beyond this point, stars begin their rapid post-AGB evolution, and they quickly move 
out of the color-magnitude space of our sample. 



4.3. Brief Summary of LPV Evolution 

LMC stars on the RGB and AGB are characterized by the presence of multiple long 
periods that show increasing length and amplitude as these stars evolve. Apart from the 
presence of the long secondary period phenomenon, which appears in approximately half of 
the stars in our sample, stars initially vary non-periodically, and only later begin to pulsate 
with periods of 20-120 days on Sequences 3 and 4. The Petersen diagram and its variants 
show that stars with their primary period on either one of these sequences often have their 
secondary period on the other sequence. Furthermore, the period changes between the first 
half and last half of the MACHO light curves suggest that the amplitudes of these pairs of 
periods are very similar for stars on the inside edges of these two sequences. Similar results 
are obtained for stars which have evolved to the luminosity at which thermal pulses begin 
on the AGB (K s rs 12); these stars are usually found on S e quenc es 1 and 2 with periods 



of 45- 500 days. This supports the arguments of IWhitelock! (119861 ) and lLattanzio fc Wood 



(120031 ). that a star excites different pulsation modes in turn as it evolves. At the points 
where the dominant modes switch, we observe pulsation in multiple periods with near equal 
strengths. The increase in the luminosity of the maximum of the luminosity function above 
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the tip of the RGB also lends support to this argument. The highest luminosity stars on 
Sequences 1 and 2 have become carbon stars. After this point, the mass- loss rate of these 
stars increases drastically and they rapidly evolve out of our sample of luminous red stars. 

This paper utilizes public domain data obtained by the MACHO Project, jointly funded 
by the U.S. Department of Energy through the University of California, Lawrence Livermore 
National Laboratory under contract No. W-7405-Eng-48, by the National Science Founda- 
tion through the Center for Particle Astrophysics of the University of California under coop- 
erative agreement AST-8809616, and by the Mount Stromlo and Siding Spring Observatory, 
part of the Australian National University. 

This work was performed under the auspices of the U.S. Department of Energy by 
Lawrence Livermore National Laboratory in part under Contract W-7405-Eng-48 and in 
part under Contract DE-AC52-07NA27344. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center /California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

The Condor Software Program (Condor) was developed by the Condor Team at the 
Computer Sciences Department of the University of Wisconsin-Madison. All rights, title, 
and interest in Condor are owned by the Condor Team. 
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Table 1. Fit Parameters of the Fourier Period-Luminosity Sequences. 



Sequence 


Band 


a 

(mag/log 10 P) 




(mag) 


Amplitude a 
(peak-to-peak) 


Optical-IR 
Phase Lag 


4 RGB b 


W 


-3.8371 


17.5432 


0.08 






J 


-2.2496 


16.7511 








H 


-2.5210 


16.2389 








K s 


-2.6689 


16.2291 






4 AGB b 


W 


-4.7744 


18.7935 


0.07 






J 


-3.8972 


18.9384 








H 


-3.9638 


18.1385 








K s 


-4.1457 


18.1713 






3 RGB b 


W 


-3.0237 


16.9699 


0.08 






J 


-2.0084 


16.7729 








H 


-2.3948 


16.4997 








K 3 


-2.4976 


16.4552 






3 AGB b 


W 


-5.0891 


20.2000 


0.17 






J 


-3.9127 


19.7273 








H 


-3.9762 


18.9294 








K s 


-4.2055 


19.0728 






2 RGB b 


W 


-2.5347 


16.4754 


0.17 






J 


-1.9506 


16.9278 








H 


-2.2132 


16.5054 








K 3 


-2.2066 


16.2846 






2 AGB b 


W 


-4.8399 


20.2378 


0.41 






J 


-2.6364 


17.8648 


0.07 


14% 




H 


-3.0089 


17.6513 


0.06 


12% 




K a 


-3.6511 


18.5875 


0.04 


14% 


1 Oxygen-stars c 












amp < 0.4 


W 


-3.9231 


19.5621 


0.22 






J 


-2.9129 


19.0857 


0.09 


24% 




H 


-3.0443 


18.4948 


0.07 


23% 




K a 


-3.2104 


18.5805 


0.06 


22% 


0.4 < amp < 1 


W 


-5.7398 


23.5574 


0.64 






J 


-3.3339 


20.1339 


0.04 


3% 




H 


-3.4560 


19.5341 


0.06 


-1% 




K a 


-3.6736 


19.7285 


0.07 


-2% 


amp > 1 


W 


-5.1763 


22.5559 


1.96 






J 


-3.2996 


20.1896 


0.31 


11% 




H 


-3.4980 


19.7971 


0.31 


13% 




K a 


-3.8292 


20.2298 


0.26 


15% 


1 Carbon-stars c 












amp < 0.4 


W 


-3.3509 


18.1327 


0.28 






J 


-1.9225 


17.1136 


0.12 


3% 




H 


-2.1524 


16.5671 


0.09 


2% 




K s 


-2.5586 


16.9922 


0.08 


16% 


0.4 < amp < 1 


W 


-4.3239 


20.4330 


0.62 






J 


-2.1333 


17.7452 


0.18 


0% 




H 


-2.5860 


17.7218 


0.13 


1% 
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Table 1 — Continued 



Sequence 


Band 


a 

(mag/ log 10 P) 


/3 

(mag) 


Amplitude a 
(pcak-to-peak) 


Optical-IR 
Phase Lag 




K a 


-3.1615 


18.5259 


0.08 


-1% 


amp > 1 


W 


-3.9039 


19.7674 


1.80 






J 


-0.4514 


13.9507 


0.58 


8% 




H 


-1.5353 


15.4351 


0.49 


12% 




K a 


-2.8910 


17.9953 


0.36 


13% 


E 


W 


-2.9765 


19.1364 


0.13 






J 


-2.4438 


19.1903 


0.05 


11% 




H 


-2.7049 


18.9063 


0.06 


6% 




K s 


-2.8198 


18.9370 


0.05 


4% 


D amp < 0.2 


W 


-2.9953 


19.7903 


0.09 






J 


-2.4171 


19.6843 


0.05 


-11% 




H 


-2.5615 


19.1820 


0.05 


-11% 




K a 


-2.7983 


19.5482 


0.05 


-12% 


D amp > 0.2 


W 


-3.3347 


20.1750 


0.48 






J 


-1.8981 


18.0679 


0.08 


-7% 




H 


-2.1999 


17.9674 


0.08 


-13% 




K s 


-2.5943 


18.7166 


0.06 


-16% 



a Amplitudes listed after W are the CLEANest amplitudes for the MACHO Blue light curve. 

b The stars above and below the tip of the Red Giant Branch (K s = 12.3) were fit separately for Sequences 2, 3, and 4. 
c The oxygen and carbon stars in Sequence 1, separated using a color cut at J — K s = 1.4, were fit individually. 



Table 2. Characteristics of the Fourier Period-Luminosity Sequences. 



Sequence Number Optical/_R" s Optical/ K s Period 

of Stars a Amplitude Ratio Phase Lag Changes' 3 

4 3,943 8% ■ • • ■ • ■ 40% do not move, 36% to D, 14% to 3 

3 4,228 9% • • • ■ • ■ 46% to D, 30% do not move, 13% to 4 

2 3,917 8% 11 14% 40% do not move, 26% to D, 22% to 1 

1 5,485 12% 4-9 -2-22% 60% do not move, 16% to 2, 15% to D 

E 888 2% 2 4% 77% do not move 

D 14,438 31% 2, 8 -12, -16% 70% do not move 

One- Year Artifact 11,215 24% 

unclassified LPVs 2,658 6% 



Number of stars with cither Po or Pi on Sequence E, Sequence D, or cither Sequence E or D. 

E 1,624 4% 

D 20,805 44% 

E or D 22,261 48% 



a The percentages are with respect to all those stars with good Red and Blue CLEANest results and a 2MASS 
K s magnitude: 46,831 stars in total. 

b From the first half of the light curve to the second half; the listed movements are just those which involve at 
least 10 percent of the stars in each sequence. 



